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Abstract. We look at the magnetic field induced weak localisation peak of graphene samples 
with different mobilities. At very low temperatures, low mobility samples exhibit a very broad 
peak as a function of the magnetic field, in contrast to higher mobility samples, where the 
weak localisation peak is very sharp. We analyze the experimental data in the context of 
the localisation length, which allows us to extract, both the localisation length and the phase 
coherence length of the samples, regardless of their mobilities. This analysis is made possible by 
the observation that the localisation length undergoes a generic weak localisation dependence 
with striking universal properties. 



1. Introduction: 

Weak localisation (WL), which is the enhanced backscattering of coherent electrons in a 
disordered media [H [2], has been observed in many materials, including in thin metals and 
in two dimensional electron gases [3j. Since WL strongly depends on phase coherence, it is 
often used to probe the phase coherence properties of a device and as a tool for extracting 
information on various scattering mechanisms. More recently, WL was discussed and observed 
in graphene 01 |6j, where the existence of two degenerate Dirac valleys leads to interesting 
new WL properties. For instance, the strength of the WL effect strongly depends on inter- 
valley scattering, and WL is expected to vanish in its absence [7]. Perturbative expressions for 
the conductivity corrections in graphene due to WL have been worked out in great detail by 
McCann et al. [5j and successfully fitted to the magnetic field dependence of the resistance at 
low temperatures [6|. The dominant feature is a peak in the resistance at zero magnetic field. 
Most of the works were done in a regime of weak disorder, where the perturbative approach by 
McCann et al. is applicable and where strong localisation is absent. 

The aim of this paper is to discuss an alternate point of view to WL, which is not restricted 
to weak disorder. This was largely driven by the desire to understand experimental observations 
of WL in a large range of graphene samples, of which some were highly disordered. At low 
temperatures, the sharpness of the WL peak in clean large scale graphene is quite striking 
[8], whereas in low mobility graphene samples, a wide negative magnetoresistance behaviour is 
observed [9j. What we show in this work, is that WL can be naturally understood within the 
framework of strong localisation and applied to the experimental case of strong disorder as well 
as weak disorder. 



2. Localisation length: 

Our results follow from the observation shown in figure 1, where the localisation length (L c ) 
exhibits a dip at zero magnetic field, which in turn leads to a peak in the resistance. L c was 
extracted by computing the resistance of a disordered graphene flake using non-equilibrium 
Green's functions. The two terminal resistance is obtained by computing the Green's function 
of the honeycomb tight binding lattice iteratively as a function of the length of the system. L c 
is then determined from the exponential increase in resistance with the length of the graphene 
flake, when the length exceeds L c . Disorder is assumed to be on-site and uncorrelated. Figure 
1 shows the dip at zero field of L c for different strengths of the disorder V, where each carbon 
site has a random potential — V/2 < V{ < V/2. V is in units of the hopping element, which is 
~ 3eV in graphene. 
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Figure 1. The magnetic field dependence 
of the extracted localisation length (in units 
of lattice sites) for different values of the 
disorder V. 
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Figure 2. The magnetic field dependence of 
the relative localisation length for the disorder 
strength V = 0.3 and for different ribbon 
widths (40 to 120 lattice sites). The line is 
fitted according to equation (1). 



L c follows the same behavior regardless of the amount of dis order o r the width of the device. 
This is nicely illustrated by rescaling the magnetic field by BVW 3 Lc and then by plotting the 
relative AL C /L C = (L C (B) — L c (0))/L c (0) as a function of field. In this case all curves overlap 
and are described by the same function within numerical accuracy. The edge plays an important 
role here because, in the weak disorder limit, we always have L c ^> where W is the width of 
the system. This is because L c scales as L c <~ W for a given disorder strength, hence no matter 
how wide the device, we will always have L c ^> W in the low disorder limit. 

In the presence of edge scattering, it was shown that the correction to the conductivity is 
given by Aa ~ (1 — 1/^1 + B 2 /B^), where B c is a characteristic field determined by the device 
width and the scattering strength [5j. Since Anderson localisation is also a measure of the return 
probability, it is not surprising that the same functional form also determines L c , hence we find 
that the expression 

AL C /L C = Lc[B l~^ M - 0.95 (l - 1/y/l + BVm) (1) 

convincingly fits the numerical data shown in figure 2. However, when the disorder becomes 
very strong as measured by L c <C W, edge scattering is no longer relevant, since in this case W 
always exceeds L c . This becomes then a two-dimensional problem, where the shortest length 
scale is given by L c . In this case expression (1) is not valid anymore and we have to use a 



two-dimensional analogue. Numerically, this is illustrated in figure 3, where we show again the 
relative L c , as a function of magnetic field, but this time for stronger disorder. When comparing 
figure 3 to figure 2, the shape has changed but not the general behavior, namely a dip at zero 
field. 
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Figure 3. The magnetic field dependence of 
the relative L c for the disorder strength V — 2 
and for different widths (40 to 120). The full 
line is obtained using equation (2). 




Figure 4. The density (or Fermi energy) 
dependence, as well as the magnetic field 
dependence of L c is plotted for V = 2 and the 
width W = 120. The Dirac point corresponds 
to E F = 0. 



The functional form of the field dependence has to be modified according to the two- 
dimensional nature of the problem and can be obtained by following, for instance, Beenakker's 
approach to WL [10] . The field correction to the conductivity can be written as Act ~ 
ln(t^/r e + 1), where ~ l/(l/r^> + 1/tb), r e is the elastic scattering time, the phase 
coherence time and tb ~ £> _1 the field induced destruction of coherent backscattering. For 
the numerical simulations, we effectively work at zero temperature, which leads to an infinite 
phase coherence time. In this case, the coherent backscattering time is not bound by but is 
bound by the localisation time, r c , where L c = y/Dr c and D is the diffusion constant. Hence, 
for the simulations, has to be replaced by r c in the expression for the field dependence. Now, 
assuming again that the relative change in L c is proportional to the coherent return probability, 
i.e., the relative change of the conductivity and further assuming that r e ~ r c , we obtain 

which nicely fits the numerical data in figure 3 (the prefactor is ~ l/ln(2) which normalizes 
expression (2) at high fields). We used t c /tb — B/B c , where B c is again the characteristic 
magnetic field. In the limit, where r e <C r c , we obtain simply AL C /L C ~ ln(l + B/B c ), which 
leads to a renormalization of B c by a factor of 2 at small fields. At low temperatures and in 
the presence of strong disorder, L c is the shortest length scale and the characteristic field is 
B c ~ Ti/eL 2 c) which defines the area of a phase coherent loop, before coherent backscattering is 
suppressed when penetrated by a flux quantum. Indeed, equation (2) fits the numerical results 
for the relative L c , when the field is rescaled by B • L 2 C) as shown in figure 3. For small disorder, 
on the other hand, where L c ^> L^, the characteristic field is instead determined by and given 
by B c ~ Ti/eL 2 ^. Hence, using equation (2) we can describe the generic behaviour of both the low 
disorder (L c ^> L^) and high disorder (L c <C L^) limits. To make the connection to localisation 



as above, it is important to distinguish between the behaviour of L c and the behaviour of the 
relative L c . Indeed, L c does not follow a generic dependence on field or density as shown in 
figure 4. It is only the relative L c which shows this generic and an almost universal behaviour 
as graphed in figure 8. 

3. Experiments: 

The next step is to apply these concepts to the experimental data. We performed experiments 
on large scale graphene as well as lithographically defined Hall bars and graphene nano-ribbons, 
in addition to large (over 100 pm) single crystal grains, of which some show strong dendritic 
structures and are dubbed graphlocons. Monolayers of graphene were grown by chemical vapor 
deposition (CVD) of hydrocarbons on 25 /xm-thick commercial Cu foils. The CVD process used 
was similar to those described in previous works [8} [TT| [T2]. 




Figure 5. The measured magneto-resistance 
of a low mobility (p ~ 100cm 2 /Vs) graphene 
Hall bar device at lOOmK. Shown is the 
relative resistance and the fit to the relative 
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The fitting parameters are B c 
-- 0.4. 



2.2T and 



c2 



o.oo 



-0.02 



-0.04 



-0.06 




Figure 6. This sample had a medium 
mobility of p ~ 700cm 2 /Vs) and the fitting 
parameters are B c = 0.08T and a = 0.06. 



The overall experimental behavior can be summarized as follows: at low temperatures, low 
mobility samples, typically below 500 cm 2 /V-s show a very wide peak in the resistance at zero 
field (as shown in figure 5). With increasing mobility, the peak becomes increasingly sharper 
(figures 6 and 7). All the samples we measured show a WL peak at zero field. The relative 
change in resistance varies between 0.5% for high mobilities to 20% for low mobilities. This is also 
expected from McCann's expression, since it yields Ap/p <~ pAF(B), where AF(B) describes 
the magnetic field dependence and p is the resistivity [5]. The main difference is that equation 
(2) leads to an increased sensitivity on the sample geometry, since for a phase coherent sample 
the resistance will depend exponentially on the length at lengths longer than the localization 
length. 

To connect the data with L c , we can simply assume that the resistance is given by R ~ e L / Lc 
where L is the length of the sample. If is smaller than the length of the sample, then L needs 
to be replaced by L^. For small changes in the relative resistance, we can then write 

AR AL C L 6 
R L c L c 



(3) 



Using equations (2) and (3) we can now fit the relative resistance as a function of the magnetic 
field for different samples. The fits are shown in figures 5 to 7 for samples with increasing 
mobilities. The fits are quite good and require only two fitting parameters, B c and a ~ 1.4L^/L C , 
which are given in the corresponding figure captions. Recalling that for large devices, the limits 
are given by B c ~ %/ emm{L 2 ^ L 2 }, which can be interpolated to B c ~ Ti(L~^ 2 + L~ 2 )/e and 
which allows us to obtain L c and L^. Plugging in the numbers we find that ~ 20nm and 
L c ~ 60nm for the sample in figure 5, whereas ~ 90nm and L c ~ 2/xm for the sample in 
figure 6, and finally ~ 0.5/xm and L c ~ 65/xm for the sample in figure 7. 




Figure 7. The relative resistance change of a 
high mobility sample (/i ~ 5000cm 2 /Vs) with 
fitting parameters B c — 0.003T and a — 0.01. 
The inset is a zoom-in of the low field region. 




Figure 8. The relative L c as a function of 
the Fermi energy and magnetic field for V = 2 
and W = 120. The absence of a systematic 
dependence on the Fermi energy outside of the 
Dirac peak is quite striking and determines a 
close to universal dependence, irrespective of 
density and disorder strength. 



4. Conclusion: 

We presented an alternate way of looking at the experimental weak localisation peak, by realizing 
that the localisation length also undergoes a weak localisation effect, in the sense that there is a 
dip at zero field, which implies a peak in resistance. This allows us to connect weak localisation 
to strong localisation in a natural way. Moreover, we can associate a characteristic field, B C1 
which determines the magnetic field dependence. B c is determined by the square of the inverse 
localization length for strong disorder and by the square of the inverse phase coherence length 
for weak disorder. Indeed, we find that the magnetic field dependence of the relative localisation 
length becomes close to universal, and that all curves can be rescaled on top of each other, by 
rescaling the field dependence with B c . We used this generic localisation length dependence to fit 
our experimental data for graphene samples of very different mobilities, which not only allowed 
us to obtain a good fit, but also to extract the localisation length and phase coherence length 
independently. This method is very robust and is applicable to high and low mobility samples. 
We expect this analysis to be very generic and applicable to a large class of two dimensional 
systems and not only to graphene. 
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